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With the purpose of controlling the steady state of a dielectric nanosphere levitated within an optical cavity, 
we study its conditional dynamics under simultaneous sideband cooling and additional time-continuous mea¬ 
surement of either the output cavity mode or the nanosphere’s position. We find that the average phonon number, 
purity and quantum squeezing of the steady-states can all be made more non-classical through the addition of 
time-continuous measurement. We predict that the continuous monitoring of the system, together with Marko¬ 
vian feedback, allows one to stabilize the dynamics for any value of the laser frequency driving the cavity. By 
considering state-of-the-art values of the experimental parameters, we prove that one can in principle obtain a 
non-classical (squeezed) steady-state with an average phonon number nph ~ 0.5. 


I. INTRODUCTION 

Bringing physical degrees of freedom to the quantum 
regime is proving so difficult that quantum control is bound 
to be a multi-branched endeavour, where techniques devel¬ 
oped on different platforms and designed for different aims 
are blended together. In the context of cooling matter to the 
quantum ground state, a primary directive of quantum control, 
various techniques have come to the fore over the last twenty 
years. Prominent among them in the case of macroscopic me¬ 
chanical systems is sideband cooling, where the targeted de¬ 
gree of freedom is driven by light on a red sideband, such that 
a beam-splitting light-matter interaction is achieved and ex¬ 
citations are drained out of the system, cooling it down. On 
the other hand, the implementation of efficient indirect quan¬ 
tum measurements is another obvious way to extract entropy 
from a quantum system. This study evaluates the combined 
performance of sideband laser cooling and continuous quan¬ 
tum measurements (also known as ‘monitoring’) on a levitat¬ 
ing nanosphere, an interesting opto-mechanical system where 
both such techniques are applicable to actual experiments. 

The research field of quantum opto-mechanics, whose goal 
is to achieve control at the quantum level of massive mechan¬ 
ical oscillators, has received increasing attention in the last 
years, both for applicative and fundamental reasons d. So 
far, the main objective pursued by both theorists and experi¬ 
mentalists is the cooling of the oscillator either to its motional 
ground state or to non-classical states, such as low-number 
Fock states or squeezed states; several protocols have been 
proposed in this respect, adapted to different physical settings 
and using different control strategies 12^171. 

Here, we will focus on a particularly promising setup, 
where the opto-mechanical system corresponds to a dielec¬ 
tric nanosphere trapped inside an optical cavity HIED- As 
the nanosphere is levitating, the coupling to the environment 
is minimised; it is possible then to neglect the thermal back¬ 
ground of phonons which is typically one of the most detri¬ 
mental sources of noise in opto-mechanical systems. Varia¬ 
tions of this basic paradigm, where the nanosphere is trapped 
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by the cavity field only |[T9j[22l, by an optical tweezer within 
the cavity field 1^ . or with the help of an electromagnetic 
trap ED, have been recently proposed. On the theoretical 
side, a detailed derivation of the master equation for the quan¬ 
tum state corresponding to the nanosphere’s motion and to the 
cavity mode may be found in ||24l|25l, and allows one to study 
the time behaviour of this opto-mechanical system, as well as 
its steady-state properties. As a matter of fact, this master 
equation paves the way to the analysis of protocols combining 
time-continuous measurement and feedback operations, and it 
will be the starting point of our study. 

We have witnessed constant progress in the understanding 
of quantum filtering, i.e. of the conditional dynamics of quan¬ 
tum systems subjected to time continuous measurements EH- 
ED. In particular, as regards systems described by continuous, 
canonical degrees of freedom, diffusive dynamics described 
by multivariate Wiener increments have been characterised in 
detail, and a general framework is available ||42l|43l. Such 
diffusive dynamics correspond to conditional evolutions due 
to the monitoring of the environment through the class of so 
called general-dyne detections ED- These quantum measure¬ 
ments amount to performing homodyne detections on the en¬ 
vironment and, possibly, additional ancillary modes which are 
coupled to the environment itself via Gaussian unitary trans¬ 
formations Ell ED. We recall that the term “homodyne” de¬ 
tection refers to the projective measurement on the eigenbases 
of canonical position and momentum operators x and p. In the 
case where the overall Hamiltonian is quadratic in the canon¬ 
ical operators, and the system is linearly coupled to the en¬ 
vironment, the conditional dynamics due to general-dyne de¬ 
tections preserves the Gaussianity of the quantum state and 
thus the whole dynamics can be equivalently described by the 
evolution of first and second statistical moments only. In this 
case, it is easy to optimise different steady-state properties, 
such as entanglement, squeezing and purity, over the parame¬ 
ters characterising the detection scheme ||46H5T]| . 

The role of conditional dynamics due to indirect measure¬ 
ments, and of subsequent quantum feedback, for steering the 
quantum state of a mechanical oscillator towards either its 
ground state or a certain non-classical state has been already 
discussed in the literature. In (521, Doherty and Jacobs derive 
the effective stochastic master equation corresponding to po- 
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sition measurements obtained through adiabatic elimination 
of a cavity mode continuously monitored through homodyne 
detection, and discuss the related feedback strategies aimed 
cooling the motion of the oscillator. The possibility to observe 
mechanical squeezing via a continuous back-action evasion 
measurement was firstly proposed in O and then invesitgated 
experimentally in El El. In (5^ . the unravelling correspond¬ 
ing to direct position measurements of an oscillator interact¬ 
ing with a non-zero temperature thermal bath is considered; 
more specifically, the results obtainable in different regimes, 
corresponding to different measurement resolutions, are dis¬ 
cussed in great detail, and the ensuing stochastic master equa¬ 
tion has been seminal in the design of protocols to engineer 
thermo-mechanical squeezing dlSl. Also, the usefulness 
of discrete and repeated measurements on a coupled qubit in 
a hybrid setup has been investigated in (551, where an effec¬ 
tive dynamics able to prepare a squeezed steady-state has been 
identified. Very recently, Hofer and Hammerer (5^ have stud¬ 
ied the effect of continuous homodyne detection on the cavity 
output combined with sideband cooling in a standard opto¬ 
mechanical setup, where the mechanical oscillator interacts 
with a non-zero temperature thermal bath. First, they consider 
a single oscillator and discuss the corresponding steady-state 
average number of phonons; then, they present more complex 
and sophisticated protocols able, for example, to create entan¬ 
glement between two distant oscillators. 

As already stated above, in this manuscript we focus on the 
case of a levitating dielectric nanosphere in an optical cav¬ 
ity, described by the master equation derived in (251 . No¬ 
tice that our treatment is distinct from the existing literature in 
that the master equation of the levitating nanosphere includes 
a photon scattering term, and the measurements of both cav¬ 
ity output, via homodyne detection, and oscillator position, 
through the light scattered by the nanosphere itself, are con¬ 
sidered simultaneously. By combining these measurements 
with sideband cooling and Markovian feedback, we address 
the possibility of both cooling the oscillator towards its ground 
state and of generating quantum mechanical squeezing, that is 
sub-vacuum fiuctuations, which is a paradigmatic signature 
of non-classicality useful for quantum metrology and preci¬ 
sion sensing E21 Finally, we also address in more detail the 
experimental setup described in 1^ . where the nanosphere is 
trapped in a high finesse optical cavity: by considering state- 
of-the-art values for the experimental parameters and for the 
measurement efficiencies, we show the possibility to vastly 
improve the performances of sideband cooling both in terms 
of steady-state average number of phonons and in terms of 
generation of squeezed quantum states. 

The manuscript is organised as follows: in Sec. |I^we dis¬ 
cuss the master equation of a levitating nanosphere, and intro¬ 
duce the notation and figures of merit that will be discussed 
in the remainder of the article. In Sec. |nl| we introduce 
the stochastic master equation describing the time-continuous 
measurements and then present the results obtainable for dif¬ 
ferent values of their measurement efficiencies. In Sec. llVIwe 
discuss the performances of these protocols for a specific ex¬ 
perimental setup, while Sec. |V|concludes the paper with some 
final remarks. 


II. LEVITATING DIELECTRIC NANOSPHERE MASTER 
EQUATION 

We will consider two quantum degrees of freedom; the 
cavity electromagnetic mode and the mechanical motion of 
a trapped nanosphere, described respectively by bosonic op¬ 
erators a and b satisfying the commutation relations [a, a^] = 
[6, 6^] = 1. We can then define the corresponding position and 
momentum quadrature operators as Xc = {apc = 
—i{a — o))/^, Xm = /\/2din&Pm = 

which can be grouped in a single vector 

^ Pct Pm) • (1) 

By considering the cavity driven by a laser at frequency ujl, 
the Hamiltonian describing the interaction between the two 
modes reads 

H = cig{aa^){bb^)^ (2) 

where g is the effective coupling constant, ujm is the mechan¬ 
ical frequency and we have already transformed the Hamilto¬ 
nian to a frame rotating at the driving laser frequency ujl, such 
that A = (jOl — ujc denotes the detuning from the cavity res¬ 
onance ujc (note that we set = 1). By considering the open 
dynamics resulting from the interaction with the environment 
(i.e. the free electromagnetic modes), one obtains the master 
equation (25l, 


= -i[H,Q]+KV[a]g + rV[b + b^]g, (3) 

where 'D[0]g = OgO^ — {O^Og + gO^O)/2. The first 
term is responsible for the unitary dynamics, the second 
one describes the usual cavity loss (with total loss rate n), 
while the third one corresponds to the recoil heating due 
to photon scattering from the oscillating nanosphere (with 
decoherence rate F). The dependence and formulas for all 
the parameters entering in Eqs. 0 and ^ can be found in 
f25\ . Specifically we want to point out that the cavity loss 
parameter n = no ^ Kd is the sum of the intrinsic loss rate no 
due to the imperfections in the cavity mirrors, plus the extra 
contribution nd due to the presence of the dielectric inside the 
cavity. We also remark that we only address the control along 
one spatial direction (dictated by the harmonic trap generated 
by the optical tweezers) along the optical cavity axis and 
we will not deal with the potential technicalities involved in 
cooling the motion along the other two decoupled directions. 

By assuming that the system is prepared in a Gaussian state 
(e.g. in a thermal state) at time t = 0, at every time the dy¬ 
namics keeps the state Gaussian (see (5^ for different reviews 
on Gaussian states). As a consequence, the whole dynamics 
can be fully described by means of the first moment vector R 
and the covariance matrix cr, whose elements are defined as 

Rj = Tr[grj] 

<^jk = T^ieirjrk + rkrj)] - 2RjRk. ( 4 ) 
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Throughout the article, we will discuss the efficiency of our 
protocols by considering the effect on the mechanical os¬ 
cillator properties. Hence, it is useful to introduce the co- 
variance (sub)matrix corresponding to the oscillator quan¬ 
tum state alone, obtained by tracing out the cavity mode, 
gim) _ JycIq]. This corresponds to 

^(m) _ f ^33 <^34 A _ 2 f ^^mPm A /^\ 

“ V <^34 0-44 y ^XmPm ^Pm ) ’ 

where (Tjk are the elements of the global covariance matrix 
defined in Eq. 0, which in fact correspond to variances and 
covariances of the quadrature operators Xm and Pm- All the 
properties we are interested in can be easily obtained from the 
matrix . We will focus on three figures of merit: quantum 
squeezing, purity and average number of phonons. 

Quantum squeezing can be quantified through the mini¬ 
mum eigenvalue of as 

^ = mineig[cr*^’^^]. (6) 


In the next sections we will show the squeezing behaviour, 
by plotting ^ in dB scale, such that negative values will cor¬ 
respond to a quadrature of the mechanical oscillator having 
sub-vacuum fluctuations, and thus to a non-classical squeezed 
state. In general this quadrature will be a certain linear com¬ 
bination of position and momentum operators, and the de¬ 
tectability and usefulness of the corresponding squeezing may 
not be straightforward. This is one of the reasons why, in Sec. 


The usefulness of a certain quantum state in quantum infor¬ 
mation and quantum communication protocols often strictly 
depends on its purity, that is on how close a state p is to a 
projector \'ip){'ip\ on a single Hilbert space vector, rather than 
a statistical mixture thereof. Such single vector states are also 
known as pure states. The purity of a quantum state is defined 
as /i = Tr[p^] and it takes is maximum value /i = 1 if and 
only if p is a pure state. For Gaussian states, the purity can 
be evaluated through the covariance matrix, and in particular 
for the single-mode mechanical oscillator we can use the for¬ 
mula /i = 1/Y^Det Notice that, as we are consider¬ 

ing single-mode states, all entropies, including the von Neu¬ 
mann entropy, are monotonic functions of the purity, which 
thus fully characterize the mixedness of the quantum state. 

Finally, as we will mainly consider steady-states of the 
mechanical oscillator having zero first moments (Tv[gXm] = 
Tr[gpm] = 0), the number of phonons can also be evaluated 
directly from the covariance matrix as 

= Trk(^;.+A)|-1 

^ tr[(7M]-2 ^ 

Here, tr['] denotes to the trace of a finite dimensional matrix. 
Notice that whenever the first moments of the oscillator are 
not equal to zero, Eq. 0 gives only a lower bound on the 
actual average number of phonons nph. 


IV 


we will also focus on the position fluctuations Ax‘^ only. 


In terms of first moment vector and of the global covari¬ 
ance matrix, the master equation for the quantum state g in 
0 entails the following time evolution: 


dR 

dt 

da 

dt 


= AR, 

= Aa aAi^ - 1 - D , 


( 8 ) 

(9) 


where the matrices A (drift matrix) and D (diffusion matrix) 
can be easily evaluated from the parameters entering in ^ 
and in the interaction Hamiltonian Q ll43]| and are reported 
in the Appendix]^ The existence of a steady-state of the dy¬ 
namics (having zero first moments) can be easily discussed by 
analyzing the eigenvalues of the drift matrix A: 


stable dynamics 4=^ Re[aj] < 0 Vj (10) 

where Re[x] denotes the real part of a complex number x, and 
aj are the eigenvalues of the drift matrix A. Such a property is 
often referred to as ‘Hurwitz stability’ in the control literature. 

A plot of the stable regions for the master equation ([^ as a 
function of the detuning A and of the coupling constant g is 
pictured in Fig. As one can see, stability is obtained only 
for red-detuning (A < 0) and, in order to reach a steady-state, 
for larger values of the coupling constant g, one needs also a 
larger value of the detuning |A|. In Fig. we plot the val- 



FIG. 1. Stable (green) and unstable (gray) regions for the mas¬ 
ter equation 0. as a function of the laser detuning A and of the 
oscillator-cavity mode coupling constant g, for k, = 2u}m and 
r = cjm/lO. Notice that, for AjijJm — 0, the drift matrix A al¬ 
ways has an eigenvalue equal to zero, and thus the system cannot be 
considered strictly stable. 

ues of the purity of the oscillator and the average number of 
phonons obtainable at steady-state as a function of the detun¬ 
ing A, where the other parameters have been given plausible 
experimental values in current setups as in ; one can see 
that for large red-detunings one can cool the oscillator to a 
state with around nph = 8 phonons at steady-state. 
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FIG. 2. Steady-state values for average number of phonons (left) 
and for the oscillator purity (right), as a function of the detuning 
A (g = ujm, 1^ — and F = ct;rn/10). Notice that we plot 
only values of detuning corresponding to the stable region and that 
minimum of phonon number and maximum of purity do not exactly 
correspond (this is due to the fact that the steady-state corresponds to 
a squeezed thermal state). 


III. TIME-CONTINUOUS HOMODYNE MEASUREMENT 
OF CAVITY MODE AND OSCILLATOR POSITION 

We now consider the conditional evolution due to general- 
dyne time-continuous measurement on both the cavity mode 
and the oscillator. The corresponding stochastic master equa¬ 
tion reads 

dg = Cgdt 1-L[ae^^]g dwi + \/r]2^ ^^[6 + b^]gdw2 

( 11 ) 

where 'H[0]g = Og -b gO^ — Tr[(0 -f 0^)g]g and dwj are 
uncorrelated Wiener increments, such that dwjdwk = dt djk- 
The term ^yr]lh^l-L[ae^^] describes the effect of continuous ho¬ 
modyne on the output cavity mode with efficiency r]i , where 
the phase (j) can be adjusted by choosing the optical phase of 
the monitored quadrature operator (e.g. 0 = 0 and 0 = 7r/2 
correspond respectively to homodyning quadratures Xc and 
Pc) 1^14^ . Analogously, the term ^/r] 2 T'H[b -f- b^ describes 
the effect of continuous monitoring of the oscillator position, 
with efficiency r ]2 ll52l[5^ . 

As for the unconditional master equation 0, the dynamics 
induced by the continuous measurement here considered does 
not change the Gaussian character of the quantum state; as a 
consequence we can translate Eq. O into equations for the 
first moment vector and covariance matrix: 

dH = AH dt -b {N — dw , (12) 

^ = Aa + aA^ — Ba + D , (13) 

dt 

where dMv = {dwi^dw 2 )^ and the matrices N,B, A and D 
can be evaluated starting from the parameters entering the 
stochastic master equation (U) ll43]| and are reported in the 
Appendix It is important to observe that the Riccati equa¬ 
tion for the covariance matrix is completely deterministic, and 
yields a steady-state that can be efficiently evaluated numer¬ 
ically. On the other hand, the first moments’ evolution is 
stochastic, i.e. it depends on the outcomes of the continuous 
measurements. As a consequence, at each time the conditional 


state is a Gaussian state whose covariances and correlations 
evolve deterministically according to the dissipative dynamics 
and the kind of measurement performed, while its first mo¬ 
ments evolves randomly in the phase-space, depending on the 
values of the photocurrents. In the following we will focus on 
these conditional steady-states only. 

Although it is possible to achieve these conditional covari¬ 
ance matrices by pure filtering, i.e. recording the measurement 
outcomes (photocurrents), in order to remain in the harmonic 
trap regime, where our treatment applies, it is useful to sup¬ 
press the drift of the first moments, due to the stochastic evolu¬ 
tion, by an active feedback operation. The role of feedback is 
indeed to use the information contained in the measurement 
outcomes in order to remove the contribution given by the 
last term in Eq. {Tg, which is proportional to the Wiener 
increment dw. This can always be done by adding a linear 
feedback term in the Hamiltonian ^ with coupling constants 
proportional to the photocurrents, i.e. 

H' = H + , (14) 


where r is the vector of quadrature operators introduced in Eq. 
0 and f (t) is an optimized vector of time-dependent coupling 
constants whose values depends linearly on the continuous- 
measurement outcomes Ea. In practice, while for the cavity 
field this corresponds simply to a linear driving, in the case 
of a mechanical oscillator it can be obtained by means of a 
combination of impulses and shifts of the trapping potential 
(for a more detailed discussion of this issue see 1521). 

The first important consequence of Eqs. in) and ( p^ re¬ 
gards the stability of the opto-mechanical system. The exis¬ 
tence of a steady-state for a continuously monitored quantum 
systems has been discussed in Ea. More specifically, it is 
proven that Eq. ( p^ has a stabilizing solution if and only if 
the pair of matrices (5, A) is detectable, namely 


5xa 7 ^ 0 V xa : Axa = Ax^ with Re[A] > 0 , (15) 


that is whenever the degrees of freedom that are not strictly 
stable under the drift matrix A contribute to the measurement 
output By. We find that, for all the choices of parameters 
we have considered in our numerical simulation, whenever 
the interaction between the two bosonic modes is on (i.e. for 
g > 0), if a continuous measurement is performed, i.e. if 
r]i > 0 or 772 > 0 , the stochastic master equation satisfies 
the stability conditions. We should remark that this stability 
condition regards the covariance matrix steady-state, while in 
principle the first moments could not go to a steady-state value 
(e.g. to zero). However, as we have just stated above and dis¬ 
cussed for example in |[52l . the information obtained from the 
measurement can be used to obtain a proper steady-state for 
the quantum system with zero first moments, as the stochastic 
drift on the latter may always be canceled by Markovian linear 
feedback. 

In the following, we will concentrate on the the steady- 
state properties of the harmonic oscilaltor. As anticipated in 
the previous section, we will analyse the number of phonons, 
the purity of the state, and at the achievable quantum squeez¬ 
ing, quantified by the minimum eigenvalue of the steady-state 
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covariance matrix. We will consider different measurement 
strategies: (i) measurement of the cavity mode only (771 > 0 
and 7 / 2 = 0 ); (ii) measurement of the oscillator position only 
( 7/1 = 0 and 7/2 > 0); (iii) simultaneous measurement of the 
cavity mode and of the oscillator ( 7/1 > 0 and 7/2 > 0). 


A. Time-continuous measurement of the cavity mode 

In this subsection we investigate the properties of the 
steady-state of the oscillator in the case where no measure¬ 
ment is performed directly on the nanosphere, while the out¬ 
put of the cavity is continuously measured. The phase (j) of the 
quadrature which is monitored through homodyne detection 
is optimized for every set of parameters and for all the figures 
of merit considered. We notice that the behaviours of these 
different optimized homodyne phases for squeezing, number 
of phonons and purity as a function of the detuning param¬ 
eter are almost identical in all the cases we investigated. In 
Fig. [^we plot the steady-state average phonon number, purity 
and squeezing for a reasonable choice of the parameters enter¬ 
ing in the master equation 0- More specifically we fix these 
parameters following the experimentally reasonable assump¬ 
tions made in Ref. ||25l, for a silica nanosphere with a radius 
of 200nm. Note that very good results are obtained for all 
the values of the detuning we are considering. This is really 
important from an experimental point of view as it strongly 
relaxes the requirement to be sideband resolved in order to 
cool the nanosphere motion, a condition that in fact is partic¬ 
ularly difficult to meet for low frequency oscillators. As one 
can notice, if we are interested in cooling the oscillator, both 
in terms of number of phonons and of purity of the quantum 
state, the optimal choice is obtained either in the red or blue 
sideband for A ^ On the other hand, in order to ob¬ 

tain the largest value of quantum squeezing, one may choose 
a value of the detuning around A = — 2.5cc^, obtaining a 
non-classical state that exhibits around 3dB of squeezing. 


B. Adding time-continuous measurement of the oscillator 
position 

As a preliminary analysis, let us consider the effect of con¬ 
tinuous monitoring of the oscillator position, whilst the cavity 
output is left unobserved ( 7/1 = 0). The results for this case 
are shown in Fig. We observe that, if we want to minimize 
the average number of phonons or maximize the purity, the 
optimal performances are obtained in the case of large detun¬ 
ing I A| ^ 1, regardless of the red or blue shift of the driving 
field. This should not come as a surprise, as a large detuning 
corresponds to decoupling the oscillator from the cavity, and 
hence directly measuring an isolated degree of freedom. By 
solving the dynamics for the decoupled mechanical oscillator 
alone (i.e., for ^ = 0), one can evaluate analytically the cor¬ 
responding steady-state. One may prove that its purity simply 
depends on the measurement efficiency, as /i = ^ 7 / 2 , which 
thus univocally characterize the entropy of the steady-state; on 
the other hand quantum squeezing and number of phonons do 



FIG. 3. Results obtained via continuous optimized homodyne of the 
output cavity mode. Top: Steady-state values for average number of 
phonons. Bottom: steady-state oscillator purity (left) and quantum 
squeezing in dB scale (right). All quantities are plotted as a func¬ 
tion of the detuning A and for different values of the cavity mode 
measurement efficiency: blue, 7 / 1=0 (notice that the blue lines are 
plotted only in the squeezing and purity plots, and only for the small 
region of values where the system is stable); purple, 771 = 0.4; yel¬ 
low, 771 = 1. The other parameters are fixed as follows: 772 = 0, 
Q — Ki — ‘2iCU'fji and F — 


depend also on the ratio between the noise parameter and the 
mechanical frequency F/cc^, and their behaviour is plotted 
in Fig. I^for different values of the measurement efficiency. 
We remind the reader that larger values of F correspond to a 
large amount of scattered light from the nanosphere. On the 
one hand, this implies a larger amount of incoherent energy 
acquired by the oscillator due to the recoil heating process; on 
the other hand, it also corresponds to a large amount of in¬ 
formation available for the continuous position measurement, 
and thus to the possibility to convert such incoherent energy 
into stead-state quantum squeezing. This explain the different 
behaviour we observe in Fig. for quantum squeezing and 
number of phonons. 

As anticipated above, the results derived here for the de¬ 
coupled oscillator almost perfectly correspond to the ones re¬ 
ported in Fig. 1^ in the case of large detuning. Note that 
the small values of the number of phonons away from res¬ 
onance are indeed due to the fact that the oscillator deco¬ 
herence rate F is relatively small with respect to (in the 
range of cc^/10). For example, for a unit efficiency measure¬ 
ment (7/2 = 1), an almost pure quantum state is obtained with 
around Tiph = 0.02 phonons. The state is also squeezed, with a 
squeezing around IdB (a value compatible with the number of 
phonons obtained). Our findings show that, if direct position 
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FIG. 4 . Results obtained via monitoring the oscillator position. Top: 
Steady-state values for the average number of phonons. Bottom: 
steady-state oscillator purity (left) and quantum squeezing in dB 
scale (right). All quantities are plotted as a function of the detun¬ 
ing A and for different values of the oscillator position measurement 
efficiency: blue, 772 = 0 . 2 ; purple, 772 = 0 . 5 ; yellow, 772 = 0 . 8 ; 
green 772 = 1 . The other parameters are fixed as follows: 771 = 0 , 
Q — (jjfjii ^ — ^^TTi and r — cc^TT^y/lO. 



FIG. 5 . Results obtained through analytical calculations for steady- 
state squeezing in dB scale (left) and number of phonons (right) 
of a decoupled mechanical oscillator (g = 0 ) subjected to time- 
continuous position measurement, as a function of the decoherence 
parameter F and for different values of the measurement efficiency 
772: from top to bottom, 772 = { 0 . 2 , 0 . 5 , 0 . 8 , 1 }. 


monitoring with high efficiency were possible, feedback cool¬ 
ing would greatly outperform sideband cooling of the oscilla¬ 
tor. However, in practice, a decoupled cavity is not likely to 
be a favourable condition to work in, as the actual efficiencies 
of position measurements through scattered light are bound 
to be severely limited by a number of practical factors (one 
among all, the geometric impossibility of probing the whole 
solid angle of scattering). Quite interestingly, at lower values 
of the measurement efficiency (i.e. 772 < 0.2) or if our aim is 


FIG. 6 . Results for simultaneous optimized continuous homodyne 
of the cavity output mode (with unit efficiency) and of the oscillator 
position Xm- Top: Steady-state values for the average number of 
phonons. Bottom: steady-state oscillator purity (left) and quantum 
squeezing in dB scale (right). All quantities are plotted as a function 
of the detuning A and for different values of the oscillator position 
measurement efficiency: blue, 772 = 0 ; purple, 772 = 0 . 5 ; yellow, 
772 = 0 . 8 ; green 772 = 1 . The other parameters are set as follows: 
771 — 1 , g — cc^TTT,, K/ — ‘ 2 iCU'fji and F — 6(7777,y^lO. 


to optimize the squeezing of the steady-state (also with larger 
values of 772 ), a combination of sideband cooling and position 
measurements still yields the best results. More specifically, 
as regards quantum squeezing, the optimal detuning is again 
around A ^ — 2.5co’^, obtaining a quantum squeezing around 
IdB and 2dB. Note that, in this section, we are not taking into 
account the fact that, in principle, varying the detuning will 
change the number of photons inside the cavity, and thus also 
the effective coupling constant g (such that larger values of 
detuning correspond to lower values of g). This effect will be 
properly taken into account in Sec.|IV] 

Finally, let us consider the simultaneous monitoring of the 
cavity output mode (optimizing the phase of the homodyne 
detection and with unit efficiency: 771 = 1 ) and of the oscilla¬ 
tor position (with different efficiencies 772 ). The results for this 
case are plotted in Fig. The properties of the steady-state 
are qualitatively similar to the ones we have just discussed, 
showing the prominent role played by the oscillator position 
measurement over the other control strategies (i.e. sideband 
cooling and cavity homodyne measurement). Nevertheless, 
one observes slightly better results with respect to the unob¬ 
served cavity scenario, both in terms of steady-state phonons 
(in particular if we do not consider the large detuning regime) 
and in terms of quantum squeezing. For red-detuning with 
A —2.5uJm, one obtains the highest value of squeezing of 
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FIG. 7 . Schematic of the experimental set-up for the quantum con¬ 
trol of a levitated dielectric nano sphere within an optical cavity. The 
nanoparticle can be sideband-cooled by sending in light that is red- 
detuned from the cavity resonance. The light that leaves the cav¬ 
ity can be continuously monitored to perform generaldyne measure¬ 
ments. The light which is scattered by the particle can be collected 
to give information upon the position of the nanosphere. 


approximately 4dB. 


IV. RESULTS FOR A NANOPSHERE LEVITATED IN A 
HIGH FINESSE OPTICAL CAVITY 


In this section, we will make specific predictions on what 
could be achieved by the continuous measurement of a 
nanosphere levitated by the field of a high finesse optical cav¬ 
ity, as depicted in Fig. [7] The position and dynamics of the 
nanosphere can be directly measured by collecting the light it 
scatters or indirectly monitored through the homodyne mon¬ 
itoring of the light that leaves the optical cavity. A detailed 
description of this setup, comprising the derivation of the for¬ 
mulas we will use in the following can be found in ITTI . 

The experimental setup we consider comprises a silica 
nanosphere of radius r = 200 nm, with mass m = 7.35 • 
10“^^Kg, and a cavity, characterized by a resonance fre¬ 
quency uJcl^Ti = 2.8 • 10^^ Hz (A = 1064 nm), length L = 13 
mm, finesse T = 400000, and cavity waist w = 60 /im . The 
corresponding value for the intrinsic cavity loss parameter is 
/^o/(27r) = cj{2TV) = 29 kHz. Following the results in 
(251, we estimate that the extra loss due to the presence of the 
nanosphere inside the cavity is of the same order, ^ 
such that the total loss parameter is fixed to k, = bS kHz. The 
cavity average photon number reads 

_ ^ Pin 1 n 

2 2/i(A-b ccc)/^^/4-b 


As usual in opto-mechanical setups, the cavity-oscillator cou¬ 
pling constant g depends on the cavity average photon number 
Tie = In particular we adopt the formula 


9 


2 


hk^A^ 

3- 

ZmUJrn 


(17) 


where k = 27r/A, 


3Vs Cr - 
2Vm + 2 


(18) 


and the volumes of the nanosphere and of the cavity read 
Vs = (4/3)7rr^ and Vm = 7 tw‘^L respectively. The remaining 
parameters, and uji denote the dielectric constant and the 
driving laser frequency. Since the nanosphere is here trapped 
by the cavity field, the mechanical frequency depends on the 
average photon number ric too, as per 


UJ 


2 

m 


2hk^A 


m 


(19) 


It is important to notice that, because of the dependence of 
uum on the average photon number ric, the behaviour of the 
opto-mechanical coupling constant is modified with respect 
to standard setups, in that g ^ nl^^. 

As far as the oscillator’s decoherence rate F is concerned, 
one can follow the results shown in (221 . One observes that 
r ^ while the corresponding ratio F /uj^ ~ 0.15 - eval¬ 
uated for our specific experimental parameters - is fixed for 
every value of the detuning A (59l . In the case of zero detun¬ 
ing (A = 0), we obtain a mechanical frequency co’^o/27r 33 

kHz and ^o/27r ^ 20 kHz. 

Like in the previous section, we will consider the steady- 
state properties as a function of the detuning A. As the me¬ 
chanical frequency uOm decreases for increasing detuning, the 
corresponding zero-point motion will increase; for this reason 
besides the quantum squeezing property (which is obtained by 
considering variances renormalized to the zero-point motion 
and thus does not take into account the effect of the variation 
of the frequency uJrn)^ we will also consider the proper fluctu¬ 
ation of the position 


(5A 


I hdxl. 


rnujm 


( 20 ) 


6x‘^ being the fiuctuations of the dimensionless position op¬ 
erator Xm- Apart from taking into account the different zero- 
point fiuctuations, 6X is arguably more interesting from an 
experimental and practical point of view as it is directly ac¬ 
cessible in experiments. 

The results are depicted in Figs. [^and|^ In the first one, we 
plot the results in the case where the efficiency of the homo¬ 
dyne measurement of the cavity output is maximum (771 = 1 ), 
the phase of the homodyne is optimized for each figure of 
merit, and different values of the oscillator’s measurement ef¬ 
ficiency r ]2 are considered. As one should expect, we observe 
better performances for increasing values of 772 . In detail, by 
focusing on the average number of phonons, for small val¬ 
ues of 772 we observe that two local minima of Upu are ob¬ 
tained, one in the red- and one in the blue-detuning regimes. 
By considering an efficiency over 50%, the optimal perfor¬ 
mances are obtained in the case of large detuning where, as 
discussed above, the oscillator’s motion decouples from the 
cavity mode, and thus its steady-state properties become en¬ 
tirely dependent on the position measurement efficiency. If we 
instead focus on quantum squeezing, in this case too we obtain 
that for all values of 772 two local minima are present in the red 
and blue sideband, and the optimal working point corresponds 
to a value of the detuning near A = — 3co’rno/2. Surprisingly, 
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FIG. 8 . Results for simultaneous optimized homodyne monitoring 
of the cavity output mode (with unit efficiency) and of the oscillator 
position Xm- Top: steady-state value for the average phonon number 
of the oscillator. Bottom left: steady-state position uncertainty 6X 
(solid lines), with the gray area corresponding to sub vacuum fiuctua- 
tions. Bottom right: steady-state quantum squeezing in dB scale. All 
quantities are plotted as a function of the detuning 5 (in units of the 
mechanical frequency uJmo, and for different values of the oscillator 
position measurement efficiency: blue, 772 = 0 ; purple, 772 = 0 . 2 ; 
yellow, 772 = 0 . 5 ; green 772 = 1 . The other experimental parameters 
are fixed as described in the main text. 


observing the behaviour of the position fluctuations dX, its 
actual minimum is always reached at resonance, while sub¬ 
vacuum fluctuations are observed only in the blue sideband, 
that is for A > 0. One important remark is needed here: 
as sub-vacuum fluctuations are not obtained for most of the 
values of the detuning A, the amount of quantum squeezing 
observed in the other plot of Fig. [^necessarily corresponds to 
different quadrature operator describing the oscillator; in gen¬ 
eral, this quadrature will correspond to a certain linear com¬ 
bination of position and momentum of the oscillator, whose 
usefulness and detectability may not be straightforward. 

Similar observations are made also regarding Fig. [^ where 
we considered plausible values for the position and homodyne 
measurement efficiencies, estimated respectively at 772 = 0 . 2 , 
based on the the set-up of Ref. ED, and at 771 = 0.5 cor¬ 
responding to the homodyne set-up efficiency in l(60l . In this 
realistic scenario, we demonstrate that for all the different val¬ 
ues of A considered one can achieve a steady-state character¬ 
ized by nph < 1, with a minimum reaching nph ~ 0.4 and a 
maximum squeezing of dB —0.5 both obtained in the red 
sideband. 





A/wmO ^/^mO 


FIG. 9 . Results for simultaneous optimized homodyne monitoring 
of the cavity output mode and of the oscillator position Xm, with real¬ 
istic values of the measurement efficiencies: 771 = 0.5 and 772 = 0 . 2 . 
Top: steady-state value for the average phonon number of the oscil¬ 
lator. Bottom left: steady-state position uncertainty 5 X (solid lines), 
with the gray area corresponding to sub vacuum fiuctuations. Bottom 
right: steady-state quantum squeezing in dB scale. All quantities are 
plotted as a function of the detuning A (in units of the mechanical 
frequency c^mo). The other experimental parameters are fixed as de¬ 
scribed in the main text. 


V. SUMMARY AND CONCLUSIONS 

Our study explicitly shows the quantum control possibili¬ 
ties offered, in realistic setups, by the combined simultaneous 
monitoring of scattered as well as coherent cavity light inter¬ 
acting with a levitating dielectric nanosphere. In particular, it 
was shown that 

• time-continuous measurements of either the cavity 
mode or the oscillator position, accompanied by Marko¬ 
vian feedback, are able to stabilize the nanosphere mo¬ 
tion for ah the values of the detuning A and of the mea¬ 
surements efficiencies. 

• The addition of time-continuous homodyne monitoring 
of the cavity output plus Markovian feedback greatly 
improves the performance that one would have obtained 
with sideband cooling only. For the realistic values of 
physical parameters considered in our study, while side¬ 
band cooling would prepare a phase-insensitive steady- 
state characterized by nph ^10 phonons on average, 
the addition of continuous homodyne measurements of 
the cavity output would prepare a squeezed steady-state, 
with nph < 2 phonons. In particular, this is true for a 
large range of detuning values, which relaxes the re- 
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quirements of sideband resolution to cool down the os¬ 
cillator. This is particularly important for low frequency 
heavy oscillators which cannot be operated in the side¬ 
band resolved regime for dispersive coupling or where 
a dissipative coupling is not available (SH. 

• In terms of optimizing the purity and the cooling of the 
oscillator, sideband cooling would cease to be useful 
if one is able to directly measure the oscillator with 
very high efficiency. Nevertheless, if we take into ac¬ 
count the state-of-the-art values for these measurement 
efficiencies in an actual experimental setup, the com¬ 
bination of the two control procedures represents still 
the best choice for experimental realizations, leading 
in principle to a quantum squeezed steady-state with 
less than one phonon on average. The performances 
are further improved if a simultaneous measurement of 
the cavity mode output is carried out. 

• If we consider a state-of-the-art experimental setup 
where the dielectric nanosphere is trapped by the field 
of a high-finesse cavity, the proposed measurement 
protocols are in principle able to prepare a quantum 
squeezed state with less than nph = 1 phonons at 
steady-state. 

Over the next few years, it will be possible to perform more 
and more exhaustive time-continuous measurements on 
the outputs of interesting micro- and meso-scopic physical 
systems 16^ . It is apparent from our findings that such a 
possibility will be one of the pathways to reduce the entropic 
content of such systems, drive them to the quantum regime, 
and ultimately achieve their full or partial quantum control. 


0 and derived following the general framework pro¬ 
vided in |^1|. 

As we stated in the manuscript, the master equation 
describing the unconditional noisy evolution of the quantum 
state of the mechanical mode and of the cavity mode, yields 
the following evolution for the corresponding first moment 
vector R and for the covariance matrix a: 


^ = Act + crA^ + D . 
at 

The drift matrix A and the diffusion matrix D read 


A = 


D = 


/ -f -A 0 0 

A -f -2g 0 

0 0 0 UJm 

V -2g 0 -UJrn 0 

^ K 0 0 0 

0/^00 
0 0 0 0 
\ 0 0 0 4r 


(Al) 

(A2) 


(A3) 


(A4) 


On the other hand, if we consider simultaneous continuous 
monitoring of the output cavity mode with homodyne detec¬ 
tion, and of the oscillator position, the corresponding stochas¬ 
tic master equation O is translated for first moment vector 
and covariance matrix in: 

dR = AR dt ^ {N — aB^) dw (A5) 

^ = Aa + aA^ — + D , (A6) 

dt 


Note added: after this work was completed, a similar ap¬ 
proach was presented in 16^ for a different experimental setup 
in which the nanosphere is trapped without the help of a opti¬ 
cal cavity. 
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Appendix A: Unconditional and conditional evolution for first 
and second moments 

In this appendix, we provide the explicit form of the matri¬ 
ces entering the unconditional and conditional evolution equa¬ 
tions for first and second moments, corresponding to the Eqs. 


where dw = {dwi^dw 2 y, A has been defined just above and 
the other matrices read 


B = 


N = 


( -\/dld cos^ (f) 

— yT^sin^cos^ 

0 

0 

/ cos^ <p 

— sin (j) cos (j) 
0 

0 


A = A + NB, 
b = D- NN'^. 


sin 0 cos 

(j) 0 

0\ 

sin^ (j) 

0 

0 

0 

0 

0 

0 

V4r?2r 0 J 



(A7) 

yTji^ sin 0 cos 

0 o\ 

y7ji77 sin^ (j) 

0 0 


0 

0 0 

5 

0 

0 0 J 




(A8) 



(A9) 



(AlO) 
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